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Abstract. We enumerate the inequivalent self-dual additive codes over 
GF(4) of blocklength n, thereby extending the sequence A090899 in The 
On-Line Encyclopedia of Integer Sequences from n = 9 to n = 12. These 
codes have a well-known interpretation as quantum codes. They can also 
be represented by graphs, where a simple graph operation generates the 
orbits of equivalent codes. We highlight the regularity and structure of 
some graphs that correspond to codes with high distance. The codes can 
also be interpreted as quadratic Boolean functions, where inequivalence 
takes on a spectral meaning. In this context we define PAR/Hiv, peak- 
to-average power ratio with respect to the {/, H, N} n transform set. We 
prove that PAR/ifjv of a Boolean function is equivalent to the the size of 
the maximum independent set over the associated orbit of graphs. Finally 
we propose a construction technique to generate Boolean functions with 
low PAR/£fjv and algebraic degree higher than 2. 



1 Self-Dual Additive Codes over GF(4) 

A quantum error-correcting code with parameters [[n, k, d]] encodes k qubits in a 
highly entangled state of n qubits such that any error affecting less than d qubits 
can be detected, and any error affecting at most qubits can be corrected. A 
quantum code of the stabilizer type corresponds to a code C C GF(4) n p. We 
denote GF(4) = {0, 1, w, w 2 }, where uj 2 = uj + 1. Conjugation in GF(4) is defined 
by x — x 2 . The trace map, tr : GF(4) i— > GF(2), is defined by tr(x) = x + x. 
The trace inner product of two vectors of length n over GF(4), u and v, is 
given by u * v — J27=i ^ r ( u iW)- Because of the structure of stabilizer codes, the 
corresponding code over GF(4), C, will be additive and satisfy u * v — for 
any two codewords it, v E C. This is equivalent to saying that the code must 
be self-orthogonal with respect to the trace inner product, i.e., C C C , where 
C 1 - = {u e GF(4)" | u * c = 0,Vc g C}. 

We will only consider codes of the special case where the dimension k = 0. 
Zero-dimensional quantum codes can be understood as highly-entangled single 
quantum states which are robust to error. These codes map to additive codes 
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over GF(4) which are self-dual 0, C = C . The number of inequivalent self- 
dual additive codes over GF(4) of blocklength n has been classified by Calder- 
bank et al. PQ for n < 5, by Hohn |H] for n < 7, by Hein et al. 4 for n < 7, and 
by Glynn et al. for n < 9. Moreover, Glynn has recently posted these results 
as sequence A090899 in The On-Line Encyclopedia of Integer Sequences jHj ■ We 
extend this sequence from n = 9 to n = 12 both for indecomposable and de- 
composable codes as shown in tabled Table shows the number of inequivalent 
indecomposable codes by distance. The distance, d, of a self-dual additive code 
over GF(4), C, is the smallest weight (i.e., number of nonzero components) of 
any nonzero codeword in C. A database of orbit representatives with information 
about orbit size, distance, and weight distribution is also available [7]. 



Table 1: Number of Inequivalent Indecomposable (i„) and (Possibly) Decomposable 
(t n ) Self-Dual Additive Codes Over GF(4) 
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1,323,363 



Table 2: Number of Indecomposable Self-Dual Additive Codes Over GF(4) by Distance 
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440 


3,132 


40,457 


1,274,068 



2 Graphs, Boolean Functions, and LC-Equivalence 

A self-dual additive code over GF(4) corresponds to a graph state @] if its gen- 
erator matrix, G, can be written as G = T + uil, where T is a symmetric 
matrix over GF(2) with zeros on the diagonal. The matrix T can be inter- 
preted as the adjacency matrix of a simple undirected graph on n vertices. It 
has been shown by Schlingcmann and Werner 0, Grassl et al. 0, Glynn US], 
and Van den Nest et al. UI] that all stabilizer states can be transformed into 
an equivalent graph state. Thus all self-dual additive codes over GF(4) can be 
represented by graphs. These codes also have another interpretation as quadratic 
Boolean functions over n variables. A quadratic function, /, can be represented 
by an adjacency matrix, T, where Tij — Tj^ = 1 if XiXj occurs in /, and Tij = 
otherwise. 
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Example 1. A self-dual additive code over GF(4) with parameters [[6,0,4]] is 
generated by the generator matrix 

(u 1 1 1 \ 

w lu 2 1 u 
ujlu 2 w 1 
1 uj uj 2 1 
1 uj 1 uj 2 
\l uj 2 uj J 

We can transform the generator matrix into the following generator matrix of 
an equivalent code corresponding to a graph state, 

(w 1 1 l\ 
Owl 1 1 
OlwllO 
lllwll 
lOllwO 
\1 1 1 u;/ 



r + Lui. 



r is the adjacency matrix of the graph shown in fig. 1(a) It can also be repre- 
sented by the quadratic Boolean function fix) = x^x?, + xqx^ + xqx§ + x,\Xi + 
X\Xz + X1X5 + x 2 x 3 + x 2 x 4 + x 3 x 4 + x 3 x 5 . 




Fig. 1: The LC Orbit of the [[6,0,4]] "Hexacode" 



Recently, Glynn et al. |5I10| has re-formulated the primitive operations that 
map equivalent self-dual additive codes over GF(4) to each other as a single, 
primitive operation on the associated graphs. This symmetry operation is re- 
ferred to as Vertex Neighbourhood Complementation (VNC). It was also discov- 
ered independently by Hein et al. 0] and by Van den Nest et al. JJ. The 
identification of this problem as a question of establishing the local unitary 
equivalence between those quantum states that can be represented as graphs 
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or Boolean functions was presented by Parker and Rijmen at SETA'01 |12| . 
Graphical representations have also been identified in the context of quantum 
codes by Schlingemann and Werner [Hj and by Grassl et al. [Hj ■ VNC is another 
name for Local Complementation (LC), referred to in the context of isotropic 
systems by Bouchet |13I14| . LC is defined as follows. 

Definition 1. Given a graph G — (V, E) and a vertex v G V. Let N v C V be the 
neighbourhood of V, i.e., the set of vertices adjacent to v. The subgraph induced 
by N v is complemented to obtain the LC image G v . 

It is easy to verify that (G v ) v = G. 

Theorem 1 (Glynn et al. [5 10J). Two graphs G and H correspond to equiv- 
alent self-dual additive codes over GF(4) iff there is a finite sequence of vertices 
v u v 2 ,...,v s , such that (((G Vl ) V2 )'") Vs =H. 

The symmetry rule can also be described in terms of quadratic Boolean 
functions. 

Definition 2. If the quadratic monomial XiXj occurs in the algebraic normal 
form of the quadratic Boolean function f , then Xi and Xj are mutual neighbours 
in the graph represented by f , as described by the n x n symmetric adjacency 
matrix, T , where T^j = J}.i = 1 if XiXj occurs in f, and T^j = otherwise. The 
quadratic Boolean functions f and f are LC equivalent if 

f'( x ) = f( x ) + ^2 x j x k (mod 2), 

j,keN a 

j<k 

where a € 1 n and N a comprises the neighbours of x a in the graph representation 
off- 

A finite number of repeated applications of the LC operation generates the orbit 
classes presented in this paper and, therefore, induces an equivalence between 
quadratic Boolean functions. We henceforth refer to this equivalence as LC- 
equivalence and the associated orbits as LC orbits. If the graph representations of 
two self-dual additive codes over GF(4) are isomorphic, they are also considered 
to be equivalent. This corresponds to a permutation of the labels of the vertices in 
the graph or the variables in the Boolean function. We only count members of an 
LC orbit up to isomorphism. As an example, fig.^shows the graph representation 
of the two only non-isomorphic members in the orbit of the [[6, 0, 4]] "Hexacode". 

A recursive algorithm, incorporating the package nauty 15 j to check for graph 
isomorphism, was used to generate the LC orbits enumerated in tabled Only 
the LC orbits of indecomposable codes (corresponding to connected graphs) were 
generated, since all decomposable codes (corresponding to unconnected graphs) 
can easily be constructed by combining indecomposable codes of shorter lengths. 

Consider, (a) self-dual additive codes over GF(4) of blocklength n, (b) pure 
quantum states of n qubits which are joint eigenvectors of a commuting set of 
operators from the Pauli Group [Q, (c) quadratic Boolean functions of n vari- 
ables, (d) undirected graphs on n vertices. Then, under a suitable interpretation, 
we consider objects (a), (b), (c), and (d) to be mathematically identical. 



Spectral Orbits and PAR of Boolean Functions w.r.t. {/, H, N} 



5 



3 Regular Graph Structures 

Although a number of constructions for self-dual additive codes over GF(4) ex- 
ist |5I16| , it appears that the underlying symmetry of their associated graphs has 
not been identified or exploited to any great extent. We highlight the regularity 
and structure of some graphs that correspond to self-dual additive codes over 
GF(4) with high distance. Of particular interest are the highly regular "nested 
clique" graphs. Fig. |21 shows a few examples of such graphs. There is an upper 
bound on the possible distance of self-dual additive codes over GF(4) Codes 
that meet this bound are called extremal. Other bounds on the distance also 
exist |1I17| . Of the codes corresponding to graphs shown in fig. El the [[6,0,4]], 
[[12,0,6]], and [[20,0,8]] codes are extremal. To find the "nested clique" graph 
representations, one may search through the appropriate LC orbits. Also note 
that all "nested clique" graphs we have identified so far have circulant adjacency 
matrices. An exhaustive search of all graphs with circulant adjacency matrices 
of up to 30 vertices has been performed. 

If d is the distance of a self-dual additive code over GF(4), then every vertex 
in the corresponding graph must have a vertex degree of at least d — 1. This 
follows from the fact that a vertex with degree 5 corresponds to a row in the 
generator matrix, and therefore a codeword, of weight 6+1. All the graphs shown 
in fig.|5]satisfy the minimum possible regular vertex degree for the given distance. 
Some extremal self-dual additive codes over GF(4) do not have any regular graph 
representation, for example the unique [[11,0,5]] and [[18,0,8]] codes. For codes 
of length above 25 and distance higher than 8 the graph structures get more 
complicated. For example, with a non-exhaustive search, we did not find a graph 
representation of a [[30, 0, 12]] code with regular vertex degree lower than 15. 



4 The {J, H, N} n Transform 

LC-equivalence between two graphs can be interpreted as an equivalence between 
the generalised Fourier spectra of the two associated Boolean functions. 

Definition 3. Let 

where i 2 = —1, be the Identity, Hadamard, and Negahadamard kernels, respec- 
tively. 

These are unitary matrices, i.e., Ji"t = HH ] = NN^ = I, where f means conju- 
gate transpose. Let / be a Boolean function on n variables and s = 2~t(— l)^( x ) 
be a vector of length 2". Let Sj, where j £ Z2™, be the jth coordinate of s. Let 
U = Uo <£> U\ <g> • • • ® U n -\ where Iff. £ {I, H, TV}, and ® is the tensor product 
(or Kronecker product) defined as 

((XqqB Oq\B ■ ■ ^ 

A ® B = a ioB anB ■ 
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Fig. 2: "Nested Clique" Graphs 
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Let S = Us for any of the 3 n valid choices of the 2™ x 2™ transform U. Then 
the set of 3™ vectors, S, is a multispectra with respect to the transform set, U, 
with 3 n 2™ spectral points. We refer to this multispectra as the spectrum with 
respect to the {/, H, N} n transform. (Using a similar terminology, the spectrum 
with respect to the {H} n transform would simply be the well-known Walsh- 
Hadamard spectrum). It can be shown that the {I, H, N} n spectrum of an LC 
orbit is invariant to within coefficient permutation. Moreover if, for a specific 
choice of U, S is flat (i.e., \S,\ = \Sj\, Vi, j), then we can write S = 
where /' is a Boolean function, h is any function from ZJJ to Zs, and u 4 = — 1. 
If the algebraic degree of h(x) is < 1, we can always eliminate h(x) by post- 
multiplication by a tensor product of matrices from T>, the set of 2 x 2 diagonal 
and anti-diagonal unitary matrices |18j , an operation that will never change the 
spectral coefficient magnitudes. Let M be the multiset of /' existing within the 
{I, H, N} n spectrum for the subcases where h(x) is of algebraic degree < 1. The 
{I, H, N} n -orbit of / is then the set of distinct members of M. In particular, if 
/ is quadratic then the {/, H, iV} n -orbit is the LC orbit |18| . 

Example 2. We look at the function f(x) = xqXi + xox%. The corresponding 
bipolar vector, ignoring the normalization factor, is 

8 = (-l)/W = (1,1, 1,-1, 1,-1,1,1) T . 

We choose the transform U = N (g> I (g> I and get the result 

S = Us = (v, v 7 , v 7 , v, v 7 , v, v, v 7 ) T , v 4 = — 1. 

We observe that \Si\ — 1, Vi, which means that S is flat and can be expressed 
as 

£ — v i(x x 1 +x X2+x 1 X2) + (6xo+6x 1 +6x2 + l) 

We observe that h(x), the terms that are not divisible by 4, are all linear or 
constant. We can therefore eliminate h{x), in this case by using the transform 

„ Ao\ Ao\ (v 7 o N 

We get the result 

JJg _ ^_^X X 1 +X X2+X 1 X2 

and thus f'(x) — xqx± + X0X2 + x±X2- The functions / and /' are in the same 
{I,H,N} n orbit, and since they are quadratic functions, the same LC orbit. 
This can be verified by applying the LC operation to the vertex corresponding 
to the variable xq in the graph representation of either function. 

5 Peak-to- Average Power Ratio w.r.t. {/, H , N} n 

Definition 4. The peak-to-average power ratio of a vector, s, with respect to 
the {I, H, N} n transform JT^ is 

PAR IHN (s) = 2™ max ISd 2 , where S = Us. 

WG{I,H,N} n 

vtez 2 » 
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If a vector, s, has a completely fiat {I,H,N} n spectrum (which is impossible) 
then PARjjjiv(s) = 1. If s = 2"i (1, 1, . . . , 1, 1) then PXR IHN (s) = 2". A 
typical vector, s, will have a PARihn(s) somewhere between these extremes. 
For quadratic functions, PAR/ hn will always be a power of 2. The PAR of s 
can be alternatively expressed in terms of the generalised nonlinearity |19| . 

7 (/) = 2f" x (2* - v/PAR/hw (i)) , 

but in this paper we use the PAR measure. Let s = 2~^ (— 1)^*), as before. 
When we talk about the PAR/ hn of / or its associated graph G, we mean 
PAR//fjv(s). It is desirable to find Boolean functions with high generalised 
nonlinearity and therefore low PARihn 120]' PAR/// at is an invariant of the 
{J, H, N} n orbit and, in particular, the LC orbit. We observe that Boolean func- 
tions from LC orbits associated with self-dual additive codes over GF(4) with 
high distance typically have low PARihn ■ This is not surprising as the distance 
of a quantum code has been shown to be equal to the recently defined Aperiodic 
Propagation Criteria distance (APC distance) |2()| of the associated quadratic 
Boolean function, and APC is derived from the aperiodic autocorrelation which 
is, in turn, the autocorrelation "dual" of the spectra with respect to {I,H,N} n . 
Table |21 shows PARihn values for every LC orbit representative for n < 12. 



Table 3: PARihn of LC Orbit Representatives 



n 








Number of orbits with specified PARihn 
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16 


32 


64 


128 


256 


512 


1024 


2048 
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1 






















2 


1 






















3 




1 




















4 




1 


1 


















5 




1 


2 
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14 
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52 


32 
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1 










9 
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156 


212 


60 


9 


1 








10 






1 


624 


1,753 


639 


103 


11 


1 






11 








3,184 


25,018 


10,500 


1,578 


163 


13 


1 




12 








12,323 


834,256 


380,722 


43,013 


3,488 


249 


16 


1 



Definition 5. Let a(G) be the independence number of a graph G, i.e., the size 
of the maximum independent set in G. Let [G] be the set of all graphs in the LC 
orbit ofG. We then define A(G) = maxnelG] a {H), i.e., the size of the maximum 
independent set over all graphs in the LC orbit of G. 

Consider as an example the Hexacode which has two non-isomorphic graphs in 
its orbit (see fig.0. It is evident that the size of the largest independent set of 
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each graph is 2, so A = 2. The values of A for all LC orbits for n < 12 clearly show 
that A and d, the distance of the associated self-dual additive code over GF(4), 
are related. LC orbits associated with codes with high distance typically have 
small values for A. Table 01 summarises this observation by giving the ranges of A 
observed for all LC orbits associated with codes of given lengths and distances. 
For instance, [[12, 0, 2]] codes exist with any value of A between 4 and 11, while 
[[12, 0, 5]] and [[12, 0, 6]] codes only exist with A = 4. 



Table 4: Range of Maximum Independent Set Size 



d 








Range 


of A for specified n 
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4 


5 


6 


7 


8 


9 


10 


11 


12 
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2,3 


3,4 


3-5 


3-6 


3-7 


4-8 


4-9 


4-10 


4-11 


3 
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3,4 


3,4 


3-5 


4-6 


4-7 


4-8 


4 










2 




3,4 


3,4 


3-5 


4-6 


4-7 


5 
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4 


6 
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Definition 6. Let A n be the minimum value of A over all LC orbits with n 
vertices. 

From table El we observe that A n = 2 for n from 3 to 6, A n — 3 for n from 7 to 
10, and A n = 4 when n is 11 or 12. 

Theorem 2. A n+ \ > A n , i.e., A n is monotonically nondecreasing when the 
number of vertices is increasing. 

Proof. Consider a graph G — (V, E) with n + 1 vertices. Select a vertex v and 
let G' be the induced subgraph on the n vertices V r \{^}. We generate the LC- 
orbit of G' . The LC operations may add or remove edges between G' and v, 
but the presence of v does not affect the LC orbit of G'. The size of the largest 
independent set in the LC orbit of G' is at least A n . This is also an independent 
set in the LC orbit of G, so A n+ i > A n . □ 

A very loose lower bound on A n can also be given. Consider a graph contain- 
ing a clique of size k. It is easy to see that an LC operation on any vertex in the 
clique will produce an independent set of size k — 1. Thus the maximum clique 
in an LC orbit, where the largest independent set has size A, can not be larger 
than A + 1. If r is the Ramsey number R(k, k + 1) [2*T]. then it is guaranteed 
that all simple undirected graphs with minimum r vertices will have either an 
independent set of size k or a clique of size k + 1. It follows that all LC orbits 
with at least r vertices must have A > k. Thus A n > k for n > r. For instance, 
i?(3, 4) = 9, so LC orbits with at least 9 vertices can not have A smaller than 3. 

For n > 12, we have computed the value of A for some graphs corresponding 
to self-dual additive codes over GF(4) with high distance. This gives us upper 



10 Lars Eirik Danielsen and Matthew G. Parker 



Table 5: Upper Bounds on A. 



n 


13 


14 


15 


16 


17 


18 


19 


20 


21 


A n < 


4 


4 


5 


5 


5 


6 


6 


6 


9 



bounds on the value of A n , as shown in table [5] The bounds on A13 and A14 are 
tight, since A12 = 4 and A„+i > A n . 

For n = 10, there is a unique LC orbit that satisfies, optimally, A = 3, 
PAR/ # n = 8 and d = 4. One of the graphs in this orbit is the graph complement 
of the "double 5-cycle" graph, shown in fig. [3J 




Fig. 3: The "Double 5-Cycle" Graph 



Theorem 3 (Parker and Rijmen [12 ). Given a graph G — (V,E) with a 
maximum independent set A C V , \A\ — a(G). Let s — (—l)^ x \ where fix) 
is the boolean function representation of G. Let U = &) ieA Hi Ii, i.e., the 

transform applying H to variables corresponding to vertices v £ A and I to all 
other variables. Then maxvt £ z 2 „ |<Sfc| 2 = 2 alyG \ where S = Us. 

Arratia et al. 22 introduced the interlace polynomial q{G, z) of a graph G. 
Aigner and van der Hoist [22} later introduced the interlace polynomial Q(G, z). 
Riera and Parker [21 showed that q(G, z) is related to the {/, H} n spectra of 
the quadratic boolean function corresponding to G, and that Q{G, z) is related 
to the {I, H, N} n spectra. 

Theorem 4 (Riera and Parker [24j ) . Let f be a quadratic boolean function 
and G its associated graph. Then PARjhn of f is equal to 2 deg Q( G ' z \ where 
deg Q(G, z) is the degree of the interlace polynomial Q{G, z). 

Theorem 5. If the maximum independent set over all graphs in the LC orbit 
[G] has size A(G), then all functions corresponding to graphs in the orbit will 
have PARjhn = 2 A(G) . 

Proof. Let us for brevity define P(G) = PAR IHN (s), where s = 2"f (-l) /(a!) , 
and f(x) is the boolean function representation of G. From theoremOlit follows 
that P(G) > 2 A (°). Choose H = (V,E) € [G] with a{H) = A(G). If \V\ = 1 
or 2, the theorem is true. We will prove the theorem for n > 2 by induction 
on |V|. We will show that P(H) < 2 al - H \ which is equivalent to saying that 
P(G) < 2 X ( G \ It follows from theorem^and the definition of Q(H, z) by Aigner 



Spectral Orbits and PAR ol Boolean Functions w.r.t. {/, H, N} 



11 



and van der Hoist [23] that P{H) = max{P(H\u), P(H u \u), P(((H u ) v ) u \u)}. 
(We recall that H u denotes the LC operation on vertex u of H .) Assume, by 
induction hypothesis, that P(H\u) = 2 X( - H ^ . Therefore, P(H\u) = 2 a( - K \") for 
some K\u £ [H\u]. Note that K\u S [H\u] implies K £ [H]. It must then be 
true that a(K\u) < a{K) < a(H), and it follows that P(H\u) < 2 a(H K Similar 
arguments hold for P(H u \u) and P(((H u ) v ) u \u), so P(H) < 2 a( - H \ □ 

As an example, the Hexacode has A = 2 and therefore PARjhn = 2 2 = 4. 

Corollary 1. Any quadratic Boolean function on n or more variables must have 
PAR IHN >2 A «. 

Definition 7. PARih is the peak-to-average power ratio with respect to the 
transform set {I,H} n , otherwise defined in the same way as PARihn- 

Definition 8. PARi is the peak-to-average power ratio with respect to the infi- 
nite transform set {U} n , consisting of matrices of the form 



where i 2 — —1, and 9 and <fi can take any real values. {U} comprises all 2 x 2 
unitary transforms to within a post-multiplication by a matrix from T>, the set 
of 2 x 2 diagonal and anti-diagonal unitary matrices. 

Theorem 6 (Parker and Rijmen [12 ). If s corresponds to a bipartite graph, 
then PARi{s) = PAR IH (s). 

It is obvious that {I,H} n C {I,H,N} n C {U} n , and therefore that PAR IH < 
PAKjh n < PAR; . We then get the following corollary of theorems El and EJ 

Corollary 2. If an LC orbit, [G], contains a bipartite graph, then all functions 
corresponding to graphs in the orbit will have PARi — 2 A ( G ) . 

Thus, all LC orbits with a bipartite member have PARjhn — PAR;. Note that 
these orbits will always have PAR; > 2 1" '1 an d that the fraction of LC orbits 
which have a bipartite member appears to decrease exponentially as the number 
of vertices increases. In the general case, PAR/ hn is only a lower bound on 
PAR;. For example, the Hexacode has PARihn — 4, but a tighter lower bound 
on PAR; is 4.486 [Ej- (This bound has later been improved to 5.103 (2SJ-) 

6 Construction for Low PAR IHN 

So far we have only considered quadratic Boolean functions which correspond 
to graphs and self-dual additive codes over GF(4). For cryptographic purposes, 
we are interested in Boolean functions of degree higher than 2. Such functions 
can be represented by hypergraphs, but they do not correspond to quantum 
stabilizer codes or self-dual additive codes over GF(4). A non-quadratic Boolean 
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function, f(x), can, however, be interpreted as a quantum state described by 
the probability distribution vector s = 2~?(— A single quantum state 
corresponds to a quantum code of dimension zero whose distance is the APC 
distance |2()j . The APC distance is the weight of the minimum weight quantum 
error operator that gives an errored state not orthogonal to the original state 
and therefore not guaranteed to be detectable. 

We are interested in finding Boolean functions of algebraic degree greater 
than 2 with low PAR/ hn, but exhaustive searching becomes infeasible with 
more than a few variables. We therefore propose a construction technique for 
nonquadratic Boolean functions with low PAR/// n using the best quadratic func- 
tions as building blocks. Before we describe our construction we must first state 
what we mean by "low PAR/ hn"- For n = 6 to n = 10 we computed PARjhn for 
samples from the space Z 2 , to determine the range of PAR/// n we can expect 
just by guessing. Table summarises these results. If we can construct Boolean 
functions with PARjhn lower than the sampled minimum, we can consider our 
construction to be somewhat successful. 



Table 6: Sampled Range of PAR/// at for n = 6 to 10 



n 


6 


7 


8 


9 


10 


# samples 


50000 


20000 


5000 


2000 


1000 


Range of PARihn 


6.5-25.0 


9.0-28.125 


12.25-28.125 


14.0625-30.25 


18.0-34.03 



Parker and Tellambura |26I27| proposed a generalisation of the Maiorana- 
McFarland construction for Boolean functions that satisfies a tight upper bound 
on PAR with respect to the {H, N} n transform (and other transform sets), this 
being a form of Golay Complementary Set construction and a generalisation of 
the construction of Rudin and Shapiro and of Davis and Jedwab 2§]. Let p(x) 
be a Boolean function on n = X^=o tj variables, where T — {to, t%, . . . , tt-i} 
is a set of positive integers and x 6 ZJ?. Let yj £ Z 2 J , < j < L, such that 
x = yo x Hi x • ■ ■ x Ul—1- Construct p(x) as follows. 

L-2 L-1 

p( x ) = ](yjhj(yj+i) + 9j(vj), (!) 

j=0 j=0 

where 0j is a permutation: Z 2 3 — > Z 2 3+1 , is a permutation: Z 2 3+1 — > Z 2 3 , and gj 
is any Boolean function on tj variables. It has been shown |57| that the function 
p(x) will have PARz/n < 2* max , where £ max is the largest integer in T. It is 
helpful to visualise this construction graphically, as in fig. In this example, 
the size of the largest partition is 3, so PAR//tv < 8, regardless of what choices 
we make for 6 j , jj , and gj . 

Observe that if we set L — 2, t = to = t\, let #0 be the identity permutation, 
and go — 0, construction (|TJ reduces to the Maiorana-McFarland construction 
over 2t variables. Construction can also be viewed as a generalisation of the 



Spectral Orbits and PAR of Boolean Functions w.r.t. {/, H, N} 



13 




"path graph", f(x) — xoXi+xix 2 + - ■ - + x n -2X n -i, which has optimal PAR with 
respect to {H, TV}". Unfortunately, the "path graph" is not a particularly good 
construction for low PAR/ hn ■ But as we have seen, graphs corresponding to self- 
dual additive codes over GF(4) with high distance do give us Boolean functions 
with low PARihn- We therefore propose the following generalised construction. 



L-1 L-1 



L-1 



(2) 



=0 j=i+l 



where J^- is either a permutation: Z 2 — > Z 2 J , or r^j — 0, and gj is any Boolean 
function on tj variables. It is evident that r can be thought of as a "generalised 
adjacency matrix", where the entries, fjj, are no longer or 1 but, instead, 
or permutations from Zo* to Z 2 3 . Construction JJJ then becomes a special case 



where rjy = except for when j 



1 (i.e., the "generalised adjacency ma- 



trix" of the "path graph"). In order to minimise PARihn we choose the form 
of the matrix r according to the adjacency matrix of a self-dual additive code 
over GF(4) with high distance. We also choose the "offset" functions, gj, to be 
Boolean functions corresponding to self-dual additive codes over GF(4) with 
high distance. Finally for the non-zero rjy entries, we choose selected permuta- 
tions, preferably nonlinear to increase the overall degree. Here are some initial 
results which demonstrate that, using (J2J, we can construct Boolean functions 
of algebraic degree greater than 2 with low PARihn- (We use an abbreviated 
ANF notation for some many-term Boolean functions, e.g. 012, 12, is short for 

X XiX 2 + X\X 2 + Xq.) 

Example 3 (n=8). Use the Hexacode graph / = 01, 02, 03, 04, 05, 12, 23, 34, 
45, 51 as a template. Let to — 3, t\ = t 2 = t 3 = t 4 = t 5 = 1. (See fig. [5J) We use 
the following matrix r. 



r = 



/o 

3 
4 
5 
6 



02,1 

4 


7 



02,1 
3 

5 





02,1 

4 

6 




02,1 


5 

7 



02, 1\ 
3 


6 

J 



Let go(yo) — 01, 02, 12 and all other gj any arbitrary affine functions. Then, 
using © to construct p(x) we get p(x) = 023, 024, 025, 026, 027, 01, 02, 12, 13, 
14, 15, 16, 17, 34, 37, 45, 56, 67. Then p(x) has PAR IHN = 9.0. 
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Fig. 5: Example ol Construction with low PAR/£rjv 



Example 4 (n = 8). Use the Hexacode graph / = 01, 02, 03, 04, 05, 12, 23, 34, 
45, 51 as a template. Let to — 3, t\ = t 2 = t 3 = i 4 = i 5 = 1. (See fig. EJ) We use 
the following matrix r. 



r = 



/0 02,1 

3 

4 4 

5 

6 
\7 7 



12,0,1,2 01,02,12,1,2 01,02,12 02,12,1, 2\ 




4 

G 




Let go(yo) — 01, 12 and all other gj any arbitrary affine functions. Then, using J5J 
to construct p(x) we get p(x) = 015, 016, 023, 025, 026, 027, 124, 125, 126, 
127, 01, 04, 12, 13, 14, 15, 17, 24, 25, 27, 34, 37, 45, 56, 67. Then p(x) has 
PAR /HA r = 9.0. 

Example 5 (n = 9). Use the triangle graph / = 01,02,12 as a template. Let 
to = t\ = t<z = 3. (See fig. El) Assign the permutations 

r ,i=r 0)2 = (12,0,1,2)(01,2)(02,1,2) 
Ji, = (34,5)(35,4,5)(45,3,4,5) 
A, 2 = (45,3,4,5)(34,5)(35,4,5) 
r 2 , = (68,7,8)(78,6,7,8)(67,8) 
/V = (78,6,7,8)(67,8)(68,7,8) 

Let g {y ) = 01,02,12, 51(2/1) = 34,35,45, and g 2 (y 2 ) = 67,68,78. Then, us- 
ing © to construct p(x) we get, p{x) = 0135, 0178, 0245, 0267, 1234, 1268, 
3467, 3568, 4578, 014, 015, 016, 017, 018, 023, 024, 025, 028, 034, 068, 125, 127, 
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128, 134, 145, 167, 168, 234, 235, 245, 267, 268, 278, 348, 357, 358, 378, 456, 457, 
458, 468, 478, 567, 568, 578, 05, 07, 08, 13, 14, 17, 23, 25, 26, 28, 36, 37, 38, 46, 
56, 58, 01, 02, 12, 34, 35, 45, 67, 68, 78. Then p(x) has PAR IHN = 10.25. 




Fig. 6: Example of Construction with low PAKihn 

The examples of our construction satisfy a low PAR/ hn ■ Further work should 
ascertain the proper choice of permutations. Finally, there is an even more ob- 
vious variation of construction J5J), suggested by the graphs of fig. where the 
functions gj are chosen either to be quadratic cliques or to be further "nested" 
versions of construction We will report on this variation in a future paper. 
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